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1 Introduction 



Heavy quark production in hadronic collisions is one of the most interesting testing 
grounds of Quantum Chromodinamics. The fusion reactions g g Q Q and q q —>■ 
Q Q are expected to be dominant in heavy quark, Q, production and, as a matter 
of fact experimental data seems to reasonably agree with perturbative QCD (pQCD) 
calculations which today are available up to Next to Leading Order (NLO) [Ij. 

Once a heavy quark is produced, it has to hadronize to produce the observed, 
hadronic, final state. On this respect, it was expected that the factorization the- 
orem be valid and hadronization proceeded through the fragmentation mechanism. 
Thus, heavy hadron production can be separated into the perturbatively calculable 
hard scattering and gluon dynamics from the non-perturbative bound state dynamics 
contained in the process independent hadron structure, expresed through the corre- 
sponding parton distribution functions (PDF), q{x, Q"^) and g{x, Q^), and the jet frag- 
mentation functions, D/i/q(2;, Q^) 0. Literally speaking, the factorization assumption 
predicts strict independence of the heavy quark hadronization from the production 
process. Thus, no flavor correlation should exist between initial and final states. 

However, there exist copiuos experimental information on charm hadron produc- 
tion which contradicts the above hypothesis [21EIIHE1E1I3EIIH]- In fact, there was 
observed an excess in the charm hadron production at large values oi xp (~ 2pi/y/s) 
and a correlation between the leading charm hadrons with the projectile quantum 
numbers. This suggests the presence of other hadronization mechanisms which must 
be relevant at large values of xp and low p^, as long as fragmentation of the charm 
quark should not produce either flavor correlations between initial and final states, 
neither charm hadrons at large values of xp- In particular, flavor correlation suggests 
hadronization mechanims by which projetile spectators produced at small p'^ recom- 
bine with charm quarks produced either perturbatively in the hard QCD process, or 
charm quarks, though perhaps of a non-perturbative nature, already present in the 
structure of the beam particles. 

From a theoretical point of view, models have been proposed to account for the 
enhancement of charm hadron production at large xp and flavor correlations. Among 
them we can mention the model of Ref [10] , in which a charm quark produced pertur- 
batively recombines with the debris of the projectile, the intrinsic charm model 
in which a particular Fock state of the projectile containing charm quarks breaks in 
the collision giving thus the desired flavor correlation between initial and final par- 
ticles, recombination type models [I2] in which charm quarks already present in the 
projectile structure recombine with light quarks, and models based in the Dual Par- 
ton Model and Dual Topological Unitarization [ISJ in which both, the heavy quark 

is the momentum fraction of partons inside the initial hadron, z is the momentum fraction of 
the heavy quark in the final hadron. is the momentum transfer. 

leading particle is defined as the one which shares valence quarks with the colliding hadrons 
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production and hadronization, are treated on a non-perturbative basis. Likewise, 
most recently two new approachs based in recombination have been presented. In 
the first, recombination of charm and light quarks, both being part of a hard QCD 
diagram, recombine to produce a charm hadron [H]. In this scheme, recombination is 
thought as a higher order correction to hadronization since hard scattering diagrams 
contributing to this process are of NLO or higher. The second approach [15] involves 
a modification of the usual recombination prescription to produce an enhancement in 
the central (low xp) production of charm hadrons. 

All the above models have been more or less successful in reproducing the main 
features of charm hadron production, with possibly the only exception of the intrinsic 
charm model [11], which seems to be excluded [16] by recent experimental data on 
charm baryon production in 7r~N interactions by the E791 [2] and SELEX [3] Col- 
laborations. However, it is important to remark that, although some models are able 
to reproduce experimental data on production asymmetries, they cannot reproduce 
simultaneously data on production asymmetries and differential cross sections. This 
is the case of the intrinsic charm model, as noted in Ref. [16]. This shows that in 
order to make a meaningful comparison among models and experimental data, both, 
the differential cross section and the production asymmetry have to be taken prop- 
erly into account. As a matter of fact, we have two out of three quantities which are 
independent, namely, the differential cross sections of both, particle and antiparticle, 
or one of the cross sections and the production asymmetry. 

In what follows we shall analyse the main features of perturbative charm pro- 
duction in hadron-hadron collisions followed by the study of possible hadronization 
mechanisms. Later, we will compare model results with available experimental data. 
The last section will be devoted to conclusions. 

2 Brief review of perturbative charm production 

In the parton model, charm quarks are produced through the interaction of partons 
in the initial hadrons. The differential cross section as a function of Xi;' is given by [17] 



where Hab is a function containing information on the structure of the colliding 
hadrons a, b, and on the hard QCD process which produces the charm quarks. At 
LO, the funtion Hab reads 



d(Jc{c) _ 1 
dxp 2 



^/s / HabiXa,Xb,IJ.F,f4i)—dppdy, 
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Figure 1: Feynman diagrams for the elementary cross sections entering in eqs. ([3]) 
and (HD at LO. 

with Xa and Xf, being the parton momentum fractions, /il^) and g{x,^'jp) the 
quark and gluon momentum distributions in the colhding particles, s the 
square of cm. energy of the partonic system and np and are the factorization 
and the renormalization scales respectively. In eq. ([1]), is the squared transverse 
momentum of the produced c-quark, y is the rapidity of the c quark and E the energy 
of the produced c-quark. The sum in eq. ([2]) runs over i = u,u, d, d, s, s. 

The elementary cross-sections for g — g anihilation and gluon fusion, da/dt \qq and 
da/dt \gg respectively, at LO are given by [T7t ITS]. 



da na"^ (fij^) cosh (Ay) + ml/ m 
'£ 9mt [1 + cosh {Ay)f 



da TTa'l Scosh (Ay) 



di '^^ 96m4 [1 + cosh {Ay)f 



, , ^ , 2mf. 2m! 
cosh {Ay) + -rf + 



(3) 



(4) 



where Ay is the rapidity gap between the produced c and c quarks and ml = ml+p^. 
The Feynman diagrams involved in the calculation of eqs. ([3]) and (jlj) are shown in 

Fig. m 

As shown in eqs. ([3}|1]), at LO the only dependence on in the elementary cross 
sections is in (fJ^'ji)- At this order, the value of the renormalization scale is fixed by 
the requirement that the propagators in diagrams of Fig. [1] be off-shell by a quantity 
of at least m^. So, it is common to use fi^ = m^. 
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Concerning the factorization scale /if, it is the scale at which the initial hadrons 
"see" one another in the collision. Then, in the above sense, fip gives the quark and 
gluon content of the initial hadrons at the time of the collision. However, differently 
to the case of the renormalization scale ^r, there is no guiding principle to help fixing 
fip- Along this work, and in order to avoid problems related with flavor excitation di- 
agrams containing heavy quarks [HIIH], we shall use fip = 1 GeV, below the threshold 
for charm production and consistent with the sum over light flavors in eq. ([2]). 

It is also known that NLO and LO calculations only differ by a /T-factor of the 
order of 2 — 3 [19], meaning that calculations can be done consistently at LO and then 
multiplied by the corresponding i^-factor. NLO calculations also show a tiny c — c 
asymmetry which is too small to produce any effect after hadronization. Note 
however that, although the small c — c asymmetry arising at NLO has the same sign 
that the observed — D~ production asymmetry (i.e. c and D~ favored over c and 
production), it has the opposite sign for baryon production. 

3 Hadronization mechanisms 

There are basically two different models for charm quark hadronization, namely, re- 
combination and fragmentation. In the first, the c (c)-quark produced in a pQCD 
process joints to the debris of the initial particles to form the final charm hadron ^10] 
while in the second, the c (c)-quark fragments to the final charm hadron leaving a 
string of quarks behind it [20]. Variation of these processes can also be found in the 
literature [21]. However, there are some features which are common to any class of 
hadronization processes: hadrons are colorless, which means that hadronizing quarks 
must be in a color singlet state and, at the end of hadronization, no free quarks could 
exist any more. These two basic requirements are the consequence of confinement in 
QCD. Consequently, whatever the hadronization process be, a color string must be 
formed among the hadronizing quarks which has to have the correct color quantum 
numbers and the exact number of quarks to ensure that no free quarks remain at the 
end of the process and colorless hadrons are formed. These lead us to the following 
clasification: i) Short color strings, responsible for recombination processes, in which 
a q — q pair or a qqq forms the final hadron without further hadron emission, and 
a) Large color strings, producing the final hadron by fragmentation being it accom- 
panied by the emission of light (mostly pion) mesons. Note that, in addition, for 
each case there exist also two different types of color strings, namely, the meson and 
the baryon-like strings, being them characterized by the fact that the first is formed 
among a quark and an antiquark and the second one between a quark and a diquark, 
as shown in Fig. [2j The first will produce predominantly mesons while the later has 
to produce at least one baryon in order to conserve the baryon number. 

In what follows, we will analyze both, fragmentation and recombination, having 
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Figure 2: Allowed color string configurations: meson- like string (a) and baryon-like 
string (b). 

in mind their applicability to heavy quark hadronization process. 
3.1 Fragmentation 

A typical quark configuration leading to fragmentation is shown in Fig. [31 Formation 
of color strings among the charm quarks and the remnants of the colliding hadrons 
would give rise to the production of open charm, i.e. D mesons and eventually charm 
baryons. 

This kind of contribution to hadronization is usually modeled by the convolution of 
the c-quark differential cross section with a Peterson Fragmentation Function, Dh/ci 



dan f dz dac(c) 



dxp 
z 



Xp 
X 

N 

Dh/c{z) = — ^, (5) 

z(l-^-j^ 

V Z 1 — 2 

where dac{c)/dx is given in eq. ([I]). Along this work we will use e = 0.022 for D meson 
production and e = 0.06 for Ac. These values for e have been extracted from charm 
hadron production in — e~ interactions [22j. Note that the fragmentation function 
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Figure 3: Typical configuration leading to production in p — p interactions by 
color string fragmentation. 

is used to model the production of a hadron containing a heavy quark at the end of 
the color string, regardless of what happens with the remaining of the color string. 



3.2 Recombination 

A typical quark configuration leading to recombination is shown in Fig. |H The 
contribution of the recombination mechanism can be estimated by means of 

da Vs f ^ da* dx^dz 

XyZ - — — n[Xv, z,xf) , (oj 



dxp 2 J dxvdz Xv z 

where R{xv, z,xf) is the recombination function, for which we shall use 

R (x^, z, xf) = P^^^6 (1 - , (7) 

Xp \ Xf / 

with (3 a normalization parameter which has to be fixed from experimental data. In 
eq. z* = 2Ecl^ and 

da* . 2\ d.a 



da _ fW /"l-^" Hai,{Xa,Xb,IJ,l,IJ,ji) 



dXydz Jo Jz+/{l-z^) Xa — Z- 



a 5 
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Figure 4: Typical configuration leading to A+ production in p — p interactions by 
recombination. 

where Hab is given in eq. ([2]), Xy is the fraction of the momentum of the hadron a 
carried by the spectator quark g^,, z = 2pz^c/ y/s, z± = ^{z* ± z), Xb = XaZ-j [xa — 2+) 
and 



= Jrnl+pl 



{2-x, 



In eq. ([8]), is the probability density function of a quark in the initial hadron. For 
baryon production, g„ has to be replaced by x^25 I^'f) = 1vi{.x^i, ix^p) qv2{.Xy2, A*f) 

and the recombination function of eq. ([7]) by [25] 

R{x,y,z,XF) = /3 — 5— () 1 • (10) 

Xp \ Xf / 

Notice that, as the quark contains information about the structure of one of the 
colliding hadrons, this introduces a flavor correlation among the initial and final 
states. 



3.3 Direct recombination 

A second possibility is that a spectator charm quark in a flavor excitation subdiagram 
recombines with the debris of one of the initial hadrons as shown in Fig. O This 
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Figure 5: Typical configuration leading to A+ production in p — p interactions by 
spectator recombination. 

kind of contributions can be calculated along the lines of the model of Refs. [Ill [23] 
in which a charm quark being part of the sea of the initial hadron recombines with 
valence and sea quarks after the collision. Notice, however, that the factorization 
scale must be consistent with the fact that charm quarks exist in the sea of the initial 
hadron, i.e. fijp > m^. The production of charm mesons is then given by [24] 

dfjrec .xp dx-^ dX2 ^ , . ^ , , , . 

— = fJ F2{Xi,X2) R2{Xi,X2,Xf) , (U) 

axp Jo xi X2 

where Xi, i = 1,2, is the momentum fraction of the i^^ quark, F2{xi,X2) is the 
two-quark distribution function in the initial hadron and R2 {xi,X2,Xf) is the two- 
quark recombination function with (3 a normalization constant which must be fixed 
by comparison to experimental data. Baryon production can be calculated just by 
replacing the F2 and R2 by the corresponding functions for three quark recombination. 
Along this work we will use 

i 

F,{xi,...,Xi) = [n;.^i 

^ = 2,3, (12) 
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with 7 = 1 for i = 2 and 7 = —0.1 for i = 3 [12] and, for the recombination function 
we shall use its simplest version given by [211 [25] 



Ri{xi, ...,Xi,Xf) = — b\\ 

X'j \ Xp 



2,3 (13) 



4 Modeling and comparison to experimental data 

Several experiments have measured the differential cross section both as a function 
of xp and ])\ for charm meson and also baryon production. A few of them have also 
measured the production asymmetry, defined by 

_ dN^/dxp - dN^^/dxp 
^""^^ ~ dNL/dxp + dN^L/dxp ' ^ ^ 

where L is for Leading and NL is for Non Leading particles. 

Along this section we shall compare the production model presented in the pre- 
ceding section with experimental data for D meson production in vr" — Nucleon 
interactions and for A^, production in vr, p — Nucleon interactions, which are, up to 
our knowledge, most of the existing data in charm particle production. 

In order to compare the theoretical model to experimental data, we will use 



rec 



da da ^'^"^ , da 

a- — +h- — , (15) 



dxp dxp dx} 

where the first term in the right hand side of eq. f|T5|) accounts for charm hadron 
production through fragmentation and its expression is given in eqs. (|5]), and the 
second term represents the charm hadron production by the recombination of the 
perturbatively generated charm quark with spectator quarks, remnants of the beam 
particles, as given in eq. IQ. 

We do not include any contribution coming from the direct recombination of the 
debris of the initial particles (see Sec. 13. 3p because it is naturally included when NLO 
contributions to heavy quark production are taken properly into account, either by 
including explicitly NLO diagrams in the perturbative part of the process, either 
through the effective K-factor. However, we recognize that as charm quarks are in 
the limit of what is understood as a heavy quark, this mechanism is equally good to 
describe the recombination part of the hadronization process (See e.g. Refs [121 [HI 

m)- 

Coefficients a and b in eq. ffT^ where fixed by fitting experimental data. In most 
of the data, the differential cross section for particle and antiparticle production were 
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fitted independently by minimizing using the model of eq. f|T5|) . In these cases we 
used 



to obtain a linear system of equations in the unknowns a and b by requiring that 
has a minimum and then, the resulting 2x2 linear system of equations was 
solved analitically. When fits involved one of the differential cross sections and the 
asymmetry, the above procedure is not longer possible since a non linear 4x4 system 
of equations is obtained. In these cases we defined a combined the sum of the 

for the differential cross section plus the one corresponding to the asymmetry and 
used MINUIT to perform the fits. 

We have to remark that particle and antiparticle xp distributions and the pro- 
duction asymmetry form a set of three non-independent measurements. Then we 
performed fits to two of these measurements and derived the third from the other 
two. 

In the following we analyse charm meson and baryon data and compare with our 
model in eq. (ITS!) . 

4.1 production 

Charged D meson production has been measured in vr^ — Nucleus interactions at 350 
GeV/c by the WA92 Experiment g], 340 GeV/c by the WA82 Experiment [5], 250 
GeV/c by the E769 Collaboration [6] and 500 GeV/c by the E791 Collaboration f7]. 
E769 ^ and E791 [9J Collaborations also measured D*^ production in 250 and 500 
GeV/c 7i~ — Nucleus interactions respectively. 

In Fig. [6] we display the results of our fit to experimental data on production 
by the WA92 Collaboration [3|. The WA92 experiment used targets of copper and 
tungsten. As can be seen in the figures, the model describes quite well the experi- 
mental data. Our fit also shows that charm quark fragmentation and recombination 
of the charm quark produced in the hard QCD process with the debris of the initial 
beam particles is enough to describe the data. The asymmetry, as shown by the 
curve in Fig. [6l has been calculated from the curves obtained from fits to particle and 
antiparticle distributions using eq. f|T^ . 

In Fig. [7] we display the results of our fit to experimental data on production 
and asymmetry by the WA82 Collaboration [5j. Data were obtained using a tt^ beam 
with W/Si and W/Cu targets. 

Once again, as evidenced in the figures, the result of the fit describes well the 
experimental data on both, production and production asymmetry. It has to be 
noted however that in fits to the WA82 data, a small contribution from the 
recombination of the hard QCD charm quark with the debris of the initial pion is 
necessary, opposite to what happens with the WA92 data. This behavior is due to 




(16) 
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Figure 6: D~ and production by the WA92 Collaboration [1]. Left and middle 
panels: (full line) our model as in eq. f|T5l) compared to experimental data. Dashed and 
point-dashed lines: contributions from recombination and fragmentation respectively. 
Right: Production asymmetry. Model (full line) vs experimental data. 

the high value of the last experimental point, not seen in the WA92 data (See 
Fig. [6]). We also shown in the figure the asymmetry obtained from fits to WA92 |1] 
and E769 [6J data. As can be seen, both curves are similar and describe well the WA82 
data on production asymmetry. Note also that a null contribution from recombination 
to production results in an asymmetry growing fast than that shown by the full 
line in Fig. [71 giving a better agreement among data and model. 

In Fig. [8] we show the data and fits to the D*(2010) production in 500 GeV/c 
TT' Nucleus interactions by the E791 Collaboration |9j. The E791 measured the D*~ + 
D*^ particle distribution function of well as the D*^ production asymmetry. 
In order to fit the data with our model of eq. 0151) . we fitted simultaneously the 
D*^ + D*^ particle distribution and the production asymmetry. From the fits, we 
extracted the individual D*~ and D*^ particle distributions. As can be seen in the 
figure, our curves agree well with the experimental data. Data on production 
asymmetry from the same experiment [7j is superimposed on the D*^ data. As 
shown in the figure, both charged D and D* asymmetries are the same within errors. 
The E769 Collaboration has also measured the xp D*^ and distributions and 
asymmetry in 250 GeV/c vr-nucleon interactions [8]. They reported that the behavior 
of the Xp particle distributions are of the form (1 — x)^'^^"'^ and (1 — x)^'^^*^'^ for 
leading and no-leading D* mesons respectively in the region 0.1 < < 0.6. From 
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Figure 7: Left: and production by the WA82 Collaboration [5]. Left and 
middle panels: (full line) our model as in eq. f|T5l) compared to experimental data. 
Dashed and point-dashed lines: contributions from recombination and fragmenta- 
tion respectively. Right: Production asymmetry. Model (full line) vs experimental 
data. Also shown te asymmetry as fitted to the WA92 data (dashed line) and the 
one extracted from measurements of the D~ and xp distributions by the E769 
Collaboration [6] (dashed-point line). 

these fits we extracted the asymmetry shown by the dashed line in Fig. [8], which 
reproduces the same behavior of our fit to the D*^ asymmetry measured by the E791 
Collaboration. This gives support to the idea that the production mechanisms of the 
pseudoscalar and vector D*"^ mesons are the same, depending only on the flavor 
content of the initial and final particles. 

A comparison of all the data and curves on production asymmetries shows 
that the asymmetry seems to be independent of the collision energy. 

4.2 and production 

In Fig. we display data on and production by the WA92 [4j experiment and 
compare to the model of eq. ( fT5|) . As evidenced in the figure, the model describes 
quite well the experimental data. Notice that in vr^ — Nucleus interactions, the 
-D° (cm) is leading, then a positive asymmetry is expected. However, u valence quarks 
of the 7r~ can anihilate easily with u valence quarks in the Nucleus, then becoming 
unavailable to recombine with the charm quark to produce a state. Consequently 
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Figure 8: Left: D*~ + D*^ production in 7r~ Nucleus interactions at 500 GeV/c. 
Data from the E791 Collaboration [9J. Full line: our model as in eq. f|T5|) . bottom: 
Z)*~/Z}*+ production asymmetry. Full line: our fit. Data from the E791 Collaboration 
{D*: full circles [9], D^: open circles [7|). Right: D*^ and D*'^ xp distributions. Full 
line: dN/dxp as given by eq. f[T^ . Dashed line: contribution from recombination, 
point-dashed line: fragmentation. 

both and are mainly produced by charm quark fragmentation at the same rate 
and the asymmetry is largely suppressed as noted in Refs. [121 126]. Nevertheless, at 
large xp there is an excess of both and which still requires a small contribution 
from recombination. 

Table [1] lists the results of the fits to experimental data on D production and 
asymmetries. As shown in Table [H the of fits are in general larger than 
desirable. This somewhat unpleasant aspect of the fits is mainly due to two reasons, 

1) the bad quality of experimental data exhibiting large error bars in most cases and, 

2) the small number of parameters of the model, which do not allow to describe the 
details of the data but the more or less gross features of them. 

4.3 production 

In Fig. [ini we display the data and fit results for Ac production and production 
asymmetry obtained in 600 GeV/c Nucleus interactions by the SELEX Collab- 
oration [3]. We display also the E791 [2] data on A^ production asymmetry in 
7r~ — Nucleus interactions at 500 GeV/c. 

Here we fitted the A+ and A~ particle distributions and obtained the asymmetry 
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Figure 9: and D production by the WA92 Collaboration [4]. Left and middle 
panels: (full line) our model as in eq. ( |T5|) compared to experimental data. Dashed 
and point line: contributions from fragmentation and recombination respectively (See 
the text). Right: Production asymmetry. Model (full line) vs experimental data. 

from our fitting functions. As can be seen in the figures, our theoretical curves agree 
well with the experimental data. 

It is interesting to note that, although both A+ (udc) and A~ (udc) are leading in 
TT" (■U(i)-Nucleon interactions, there is a significant difference in their xp distributions. 
In fact, as noted in Ref. [16], valence u quarks in the vr" can annihilate easily with 
(valence) u quarks in nucleons but, valence d quarks of the pion cannot, thus favoring 
the A^ production over the A~ one. This is the case of diagrams c) and d) in 
Fig. [TTJ However, not only the recombination mechanism contributes to the A^ 
production asymmetry in vr" Nucleon interactions. Color string fragmentation also 
gives a contribution. Looking at the diagrams a) and b) in Fig. [11] we note that also 
in this case A~ production is disfavored. An analogue diagram to that of Fig. [TT] a), 
with the valence d quark in the vr" annihilating with a d quark in the proton and a c 
quark produced in the forward (x^ > 0) region giving rise to a A~ is suppressed by 
at least two reasons: 1) u-,^ — Up anihilation is favored over d.,^ — dp because in the 
first case both are valence quarks and 2) a c quark produced in the forward region 
has to form a baryonic color string with a ud sea diquark in the proton in order to 
produce a A~. The same situation is present in the case of Fig. [TTl d). 

Finally, in Fig. [12] we show the data and fits to A^ production and production 
asymmetry in p interactions at 600 GeV/c [3]. In order to get the theoretical 
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Experiment 


Particle 


a 


b 


xVd.o.f. 


WA82 


D- 


14269.27 ± 1116.82 


918± 140.57 


3.570 


WA82 


D+ 


18649.11 ± 1828.77 


311.99 ± 128.22 


2.368 


WA92 


D- 


292.64 ± 14.59 


9.59 ± 2.43 


3.075 


WA92 


D+ 


302.46 ± 13.92 


0. ± 1.34 


1.698 


E791 


D*- 


41.54 ± 2.14 


4.86 ± 0.54 


4.613 


E791 


D*+ 


48.40 ± 2.17 


1.97 ± 0.44 


4.613 


WA92 


DO 


630.74 ± 29.32 


4.17 ± 3.26 


3.137 


WA92 




733.51 ± 37.25 


0. ± 3.90 


4.370 



Table 1: Coefficients obtained in fits to experimental data on D meson production. 
The x^/d.o.f. of the E791 data on D* production is for a simultaneous fit to the 
D*~ + D*^ xp distribution and asymmetry. 



curves, we fitted both the A+ and A~ distributions as a function oixp and calculated 
the asymmetry from the curves obtained in fits. 

For the A~ we only use the ffist term in eq. f[T^ . neglecting any contribution from 
recombination. The reason is that in order to get a A~ baryon from recombination, 
the perturbatively produced c quark has to recombine with a ud diquark formed from 
sea quarks in the initial proton, which is less favorable than the recombination of a c 
quark with a ud valence diquark (see Fig. [T3|) . 

As in the case of A^ production in vr" — Nucleon interactions, A+ production is 
favored over the A~, having not only a different shape in their xp distributions but 
also a different global normalization, as shown in Fig. [131 This is a distinctive feature 
of charm baryon production not seen in charm meson production, where both, particle 
and antiparticle xp distributions begin at more or less the same point at ~ 0. 

Concerning the data on the production asymmetry in p Nucleus interactions 
by SELEX, it is hard to understand, however, how they obtain three points with 
asymmetry bigger than one. Note that from the definition of the asymmetry follows 
that it is bounded to be between 1 and —1. 

Table [2] shows the results of fits to experimental data on A^ production. Once 
again, as in the case of meson production, of fits are far away from their desirable 
values. However, note that in the case of baryon production, experimental data 
exhibits larger error bars than in the case of mesons. 

5 Conclusions 

We have shown that using a simple model based on perturbative QCD for charm 
quark production and on well known hadronization mechanisms, namely fragmenta- 
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Figure 10: (left) and A~ (middle) production in vr" Nucleus interactions at 600 
GeV/c. Data from the SELEX Collaboration [3j. Full line: our model as in eq. ([T5|l . 
Right: A+/A~ production asymmetry. Data from Ref. [3] (full circles) and Ref. [2] 
(open circles). Full line: result of our fit. 

tion and recombination, the available experimental data on charm hadron production 
can be well described. The idea of describing charm hadron production by means 
of fragmentation and recombination of perturbatively produced charm quarks is not 
new. It has been discussed by the first time in Ref. pilj, but never has been extensively 
tested against most of the available experimental data. 

Our analysis shows that intrinsic charm is unnecessary to describe consistently 
the experimental data on both charm particle xp distributions and charm particle 
production asymmetries. This was evident once SELEX data on A^ baryons and 
antibaryons became available, as noted in Ref. [16j. 

Parameters in the model of eq. (1151) represent the unkowns associated to the rela- 
tive fractions of the fragmentation and recombination contributions, but also include 
the uncertainties coming from the non-perturbative contributions to the hadroniza- 
tion process. It is conceivable that with more abundant and precise experimental 
data, the behavior of this parameters with respect to the reaction and reaction en- 
ergy can be fixed, allowing in this way a more predictive power of the model. 

In Refs. [121 [IS] another model including the recombination of charm quarks with 
the debris of the initial particles was considered. However, in this model, recombining 
charm quarks were considered as part of the structure of the initial particles. The 
fact that both, the model discused here and the model of Refs. [121 [IB] , are able to 
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Figure 11: Some typical diagrams for A+ and production in tt^ Nucleus, a) and 
b) A+ production by color string fragmentation. The c quark forms a baryon-like 
string with a diquark in the proton. The c quark forms a meson-like string with the 
remaining valence quark of the tt". c) and d) A^ production by recombination. 

describe well the experimental data is due to the fact that the momentum distribution 
of perturbatively produced charm quarks is similar to the sea charm quark distribution 
in hadrons. It is easy to see that this is a consequence of the fact that in both cases 
the origin of the charm quarks is gluon splitting, then its momentum distributions 
must be similar. 

In addition, we have shown that factorization is broken as long as the structure 
of the initial colliding particles has to be taken into account in order to describe the 
hadronization of charm quarks. However, it does not mean that new unknowns are 
added to the problem. There still exist a consistent way to calculate hadronization 
within the framework of the recombination model. 

Finally, wc would like to stress that not only recombination, but also fragmenta- 
tion, contributes to the observed baryon production asymetries. In recombination, 
the asymmetry is a consequence of the sharing of partons among the initial and final 
particles. In fragmentation, the asymmetry is due to the parton content of the initial 
particles. Note that, as long as initial particles are baryons and mesons, it is easier 
to have mesonic and baryonic color string, but no "anti-baryonic" ones, thus baryon 
production is favored over antibaryon production. This effect does not appear in 
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Figure 12: A+ (left) and A~ (middle) production in p Nucleus interactions at 600 
GeV/c beam energy. Data from the SELEX Collaboration [3]. Full line: our model 
as in eq. (|T5l) . Right: A+/A~ production asymmetry. Data from Ref. [3]. Full line: 
result of our fit. 

meson-antimeson production since they are formed mainly from mesonic strings. If 
antibaryons were used as initial particles, fragmentation would favor the production 
of antibaryons, producing the opposite effect. It is interesting to note also that re- 
combination gives a small contribution only noticeable at large xp (of the order of 
~ 0.6 — 0.7) to the D meson xp distributions while for baryons its contribution 
becomes important for xp ~ 0.4. 
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Figure 13: Typical diagrams for and production in p Nucleus, a) and b) 
production by recombination, c) A+ production by color string fragmentation. 



Experiment 


Particle 


a 


b 


xVd.o.f. 


71- - N 


A+ 


76511.81 ± 18661.76 


990.78 ± 303.93 


1.602 


TT- -TV 


A- 


22625.39 ± 18946.02 


1682.04 ± 478.72 


0.786 


p -N 




77904.03 ± 16361.06 


5226.22 ± 608.73 


4.077 


p -N 




32204.85 ± 19726.59 


0.00 ± 688.32 


0.093 



Table 2: Coefficients obtained from fits to data on A^ production from the SELEX 
Collaboration. 
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